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ABSTRACT 


This study investigated the effects of viewing a primary 
flight display at different retinal eccentricities on human 
manual control behavior and performance. Ten partici- 
pants performed a pitch tracking task while looking at 
a simplified primary flight display at different horizontal 
and vertical retinal eccentricities, and with two different 
controlled dynamics. Tracking performance declined at 
higher eccentricity angles and participants behaved more 
nonlinearly. The visual error rate gain increased with ec- 
centricity for single-integrator-like controlled dynamics, 
but decreased for double-integrator-like dynamics. Parti- 
cipants’ visual time delay was up to 100 ms higher at the 
highest horizontal eccentricity compared to foveal view- 
ing. Overall, vertical eccentricity had a larger impact 
than horizontal eccentricity on most of the human man- 
ual control parameters and performance. Results might 
be useful in the design of displays and procedures that in- 
crease manual control performance in critical flight con- 
ditions such as an aerodynamic stall. 


INTRODUCTION 


When manually controling an aircraft, pilots often scan a 
large visual area of the cockpit, while at the same time 
performing a manual control task. In these situations, it is 
very likely that pilots will observe the aircraft’s attitude on 
the primary flight display with their peripheral vision. Hu- 
man manual control in active compensatory tracking tasks 
is relatively well understood, and it is known that human 
controllers are able to achieve stable closed-loop perfor- 
mance by adjusting the weights put on the position and 
rate of the perceived visual cues [1]. Previous experiments 
using passive observers also showed that the perception of 
length and velocity is different in peripheral compared to 
foveal vision [2-5]. With this in mind, the aim of this re- 
search was to use a cybernetic approach to investigate the 
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effects of viewing task variables from a display peripher- 
ally at different retinal eccentricities in an active manual 
control task, by modeling human manual control behavior. 
This approach provides insights into how manual control 
performance, and the use of visual position and velocity 
information changes when peripheral vision is used in ac- 
tive manual control tasks. 


The paper adds to the literature in two ways. First, it 
investigates how human manual control parameters are af- 
fected by retinal eccentricity using a cybernetic approach. 
Second, an active control task was used, expanding on re- 
sults found in previous fundamental vision science experi- 
ments with passive observers and visual stimuli presented 
in a controlled, predefined manner. 


PREVIOUS RESEARCH 


Past studies have investigated the effects of peripheral vi- 
sion on position and velocity perception. Tynan et al. 
found that perceived velocity decreased when eccentricity 
increased. This effect reduced with higher velocities [6]. 
Moreover, both the minimum velocity perception thresh- 
old and the reaction time to motion onset increased with 
increased eccentricity. A study that tried to understand 
the neuro-physical mechanisms behind speed encoding in 
the periphery found that perceived velocity decreased with 
eccentricity only at high luminance levels, and sometimes 
the effect was inverted at very low levels of luminance [2]. 
Another study also looked at the thresholds for accelera- 
tion and deceleration detection of Gabor stimuli at differ- 
ent retinal eccentricities. Traschtz et al. found that at low 
eccentricities, humans were better at perceiving acceler- 
ation, whereas the deceleration thresholds became lower 
than those for acceleration at higher eccentricities [3]. 
Stone et al. found that speed perception is also dependent 
on the contrast of the presented stimuli [4]. 
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Thompson et al. examined the perceived length of a 
single line in peripheral view. His findings showed that 
for horizontal lines, the perceived length decreased with 
both horizontal and vertical eccentricity. For vertical lines, 
however, the perceived length appeared to decrease only 
with vertical eccentricity [5]. 

Some studies investigated the effects of peripheral vi- 
sual cues on manual control behavior and performance in 
active control tasks [7,8]. However, in these studies, pe- 
ripheral visual stimuli were always present in addition to a 
central visual display, and the effects of retinal eccentric- 
ity were not investigated. In both compensatory target- 
following and disturbance-rejection manual control tasks, 
peripheral visual cues were found to increase tracking per- 
formance, with only very limited changes in manual con- 
trol behavior. The effects of peripheral visual cues were 
more significant in target-following than in disturbance- 
rejection tasks. 

A few observations can be made after analyzing the re- 
search presented above. First, in general, length, velocity 
and acceleration determination is not as good at higher 
retinal eccentricities. Second, it is clear that speed per- 
ception is not only affected by changes in eccentricity, 
but depends on far more parameters: contrast, luminance, 
orientation, etc. Furthermore, most of the experiments 
above took place in very controlled environments, where 
the visual stimulus speed was maintained constant during 
runs, with human subjects as passive observers. In an ac- 
tive manual tracking task, the velocity of visual stimuli is 
highly variable, accelerating and decelerating depending 
on the human operator’s control inputs and the controlled 
dynamics. In addition, in such an active task, human oper- 
ators’ cognitive load is higher. Because of this, the results 
found in these passive task experiments might not directly 
apply to active manual control tasks. 


MANUAL CONTROL TASK 


A diagram showing the different components of the 
closed-loop manual control task used in this study is 
shown in Fig. 1. Here, the human operator has to activelly 
minimize the error e presented on a compensatory display, 
resembling a primary flight display (PFD), by providing 
continuous control inputs u with a joystick. These control 
inputs are transformed into aircraft pitch attitudes 6 by the 
controlled aircraft dynamics H,(s). The error e is the dif- 
ference between the actual aircraft pitch attitude 6 and a 
target pitch signal f;. Human manual control in the com- 
pensatory tracking task in Fig. 1 is typically modeled with 
a linear transfer function H,,(s), and a remnant signal n 
that captures human nonlinear behavior and noise in the 
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control loop [1]. The spectrum of the human remnant n 
has the shape of a first order low-pass filter according to 
Levison [9]. The remainder of this section provides more 
details on the different components of the control task de- 
picted in Fig. 1. 


Controlled Dynamics 


The following transfer function was used to simulate the 
controlled aircraft pitch dynamics: 


Ka 

ee s(s + wa) y) 
with kK, the gain of the controlled dynamics, and wg the 
break frequency. For values of the break frequency wa 
substantially above the crossover frequency of the human- 
operator/controlled-dynamics open loop, the dynamics are 
single-integrator-like and are perceived as easy to control. 
This is similar to controlling pitch rate. Small values of 
Wa, below the crossover frequency, represent more dif- 
ficult double integrator-like dynamics, corresponding to 
pitch-acceleration control. 


Human Operator Model 


McRuer’s crossover theorem states a human operator ad- 
justs his/her equalization dynamics to controlled dynam- 
ics such that the open-loop human-operator/controlled- 
dynamics transfer function H,H, has the form of a sin- 
gle integrator around the crossover frequency [1]. With 
the controlled dynamics of Eq. (1), human operators need 
to provide lead equalization at higher frequencies. Tak- 
ing this into account, the human operator in Fig. 1 can be 
characterized using the following transfer function: 


equalization limitations 


oor nr _— eee 


w2 
H,(s) = Kp[1 + Trs] 77" 


= 9, 
we + 2¢nWns + 8? (2) 


Equalization parameters Kk, and T’, are the human op- 
erator visual gain and lead time constant, respectively. Kk, 
is the relative weight the human operator puts on the er- 
ror signal (e), and KT’, the relative weight on error rate 
(€) in order to achieve stable control in the closed loop. 
Parameter 7,, represents the time delay associated with vi- 
sual perception, processing, and neural activation. Param- 
eters w,, and ¢, represent the neuromuscular frequency 
and damping ratio of the combined human arm/hand and 
control inceptor. 

Considering the controlled dynamics of Eq. (1), a hu- 
man operator needs to generate more lead (higher 7’, ) for 
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Figure 1 Closed-loop compensatory tracking task. 


double-integrator-like dynamics (lower values of wa) in 
order to achieve stable control. The opposite is true for 
single-integrator-like dynamics, for which less lead gen- 
eration is required. 


Target Signal 


To facilitate the identification of the linear transfer func- 
tion H,, from experimental data, the forcing function f; in 
Fig. | is typically modeled as a sum of sines, with each 
sine having a different frequency: 


Nj 
filt) = D7 Ap(k)sin[wy(k)t+ b7(k)] GB) 
k=1 


with Ap(k), w(K), and @,(k) the amplitude, frequency 
and phase of the k” sine in f;, respectively. N f Tepresents 
the number of sine waves, which was 10 in the current 
study. A summary of all forcing function properties can 
be found in Table 1. 


Table 1 Forcing function properties. 


k nf wy,tad/s Ay, deg be, deg 
1 6 0.460 6.2472 -84.774 
2 13 0.997 4.3688 -4.269 
3 27 2.070 1.9712 40.141 
4 41 3.144 1.0616 —_ -112.088 
5 53 4.065 0.7128 — -161.179 
6 73 5.599 0.4416 120.470 
7 103 7.900 0.2808 = -149.989 
8 139 10.661 0.2048 129.202 
9 174 13.346 0.1712 -38.612 

10 229 17.564 0.1456 11.127 


The sinusoid frequencies were all integer multiples n + 
of the measurement time base frequency, Ww, = 27/Tm = 
0.0767 rad/s. Ty, = 81.92 s was the measurement time 
used for the experiment. The selected integer multiples 
were used in a previous experiment and ensured that the 
ten sinusoid frequencies covered the frequency range of 
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human control at regular intervals on a logarithmic scale. 
More details on the used forcing function can be found 
in [10]. 


EXPERIMENTAL SETUP 


Method 


Apparatus 


The experiment setup is shown in Fig. 2. Participants were 
seated in front of a display and were instructed to perform 
the manual control task using a BG Systems joystick lo- 
cated on the right side. A head and chin rest were used 
to ensure that the distance from the display was identical 
for all participants. The height of the table on which the 
display and head rest were mounted could be adjusted, in 
order to accommodate different participants’ heights. The 
distance from the head rest to the display was set to 14 
inches to allow visual stimuli to be presented at predeter- 
mined retinal eccentricities. This distance allowed for the 
largest horizontal eccentricity to be approximately 30 deg, 
given the size of the monitor. 

A simplified PFD of square dimensions was presented 
at different retinal eccentricities on a 27” Apple moni- 
tor having a resolution of 2560 x 1440 pixels. The PFD 
spanned a field of view of approximately +-5 deg. On 
the PFD, a line representing the horizon divided the blue 
color at the top and the dark-brown color at the bottom. 
An aircraft symbol was fixed in the middle of the PFD. 
The error e (Fig. 1) was the difference between the hor- 
izontal lines of the fixed aircraft symbol and the moving 
horizon. A small red cross in the middle of the monitor 
was the fixation point where subjects were instructed to 
fixate at all times. The position of the PFD changed be- 
tween the different experimental conditions. The fixation 
point overlapped with the center black square of the fixed 
aircraft symbol in the foveal condition. The background 
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Figure 2 Experimental setup. 


Table 2 Experimental Conditions. 


cond. controlled PED position (POS) factor(level) 
dynamics (DYN) 
Gl velocity (SI foveal (F) DYN(1),POS(1) 
C2 velocity (SI +15 deg horizontally (R) DYN(1),POS(2) 
C3 velocity (SI) +30 deg horizontally (RR) DYN(1),POS(3) 
C4 velocity (SI -15 deg vertically (B) DYN(1),POS(4) 
C5 acceleration (DI) foveal (F) DYN(2),POS(1) 
C6 acceleration (DI) +15 deg horizontally (R) DYN(2),POS(2) 


C7 acceleration (DI) +30 deg horizontally (RR) DYN(2),POS(3) 
C8 acceleration (DI) -15 deg horizontally (B) DYN(2),POS(4) 


color of the monitor was a shade of dark grey, chosen such 
that no after-image effects would occur during the experi- 
ment. Fig. 3 shows the display with the PFD at all possible 
eccentricities (Table 2). 


Conditions 


The experiment had two independent variables: eccentric- 
ity angle of the PFD with respect to the fixation point, 
and the type of controlled dynamics. Four eccentricity an- 
gles and two types of controlled dynamics were tested in 
a full-factorial design. A summary of the total of eight 
conditions is given in Table 2. 

The four positions of the PFD are indicated by F, R, RR, 
and B in Fig. 3. F represents the foveal condition, R the 
+15 deg horizontal eccentricity condition (right), RR the 
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primary flight display 


fixation point 


Figure 3 Experimental display. 


+30 deg horizontal condition (far right), and B the -15 deg 
vertical eccentricity condition (bottom). Foveally viewing 
the control task was the baseline condition. To investi- 
gate possible trends in manual control parameters induced 
by retinal eccentricity, two different horizontal eccentric- 
ities were tested. Thirty degrees was the largest achiev- 
able eccentricity angle with the experimental setup used. 
The second horizontal eccentricity was chosen to be half 
of that. In addition, to investigate the difference between 
horizontal and vertical eccentricity, a vertical eccentricity 
was tested with the same angle as the middle horizontal 
eccentricity (i.e., 15 deg). 

The two controlled dynamics are indicated by SI and 
DI (Table 2). SI represents easier single-integrator-like 
dynamics (Kg = 10.5, wq = 6.0 s), whereas DI rep- 
resents more unstable, double-integrator-like dynamics 
(Kq = 4.9, wa = 1.08), see Eq. (1). These two controlled 
dynamics were chosen to provide insights into whether 
changes in manual control behavior induced by retinal ec- 
centricity depend on the controlled element. The break 
frequencies (wz) of the two controlled dynamics were cho- 
sen to achieve the desired variation in control difficulty, 
after which the gains (/‘q) were tuned to create a simi- 
lar control authority for both. The controlled dynamics 
were not chosen to be pure single and double integrators 
in order to estimate the lead time constant more accu- 
rately [1,9]. Furthermore, single-integrator-like dynam- 
ics are representative for an aircraft flying under normal 
conditions, whereas double-integrator-like dynamics oc- 
cur during unstable flight, in situations such as an aerody- 
namic stall. 


Participants and Procedures 


Ten participants between the ages of 22 and 58 partici- 
pated in the experiment. Four had considerable experience 
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with manual tracking tasks, one was a commercial pilot, 
four were general aviation pilots, and one was a gradu- 
ate student with no prior experience with manual tracking 
tasks. 

Prior to the experiment, participants received a briefing, 
explaining the task and how to operate the joystick. Parti- 
cipants were instructed to continuously minimize the error 
on the PFD (i.e., keep the aircraft symbol on the horizon), 
while fixating on the red cross in the middle of the display 
and giving smooth, continuous inputs. 

Each participant performed a total of 56 runs, 28 for 
each controlled dynamics. Each run lasted 90 seconds. 
To minimize the effects of adaptation from one controlled 
dynamics to another, all runs for a particular controlled 
dynamics were presented in one segment, followed by all 
runs for the other dynamics. Half of the participants per- 
formed 28 runs of conditions C1-C4 presented according 
to Latin-square design first, and then 28 runs of conditions 
C5-C8, also presented according to a Latin-square design. 
For the other half of the participants, the order of the dy- 
namics was reversed. The first eights runs for each con- 
trolled dynamics were used as training and familiarization 
runs, and were not used for data analysis. 


Dependent Variables 


The goal of the experiment was to investigate the effects 
of viewing task variables at different degrees of retinal 
eccentricity on manual control behavior. Therefore, hu- 
man control behavior parameters and performance were 
the variables of interest. 

The root mean square (RMS) of the error signal RIZS, 
was used as a measure for tracking performance. A lower 
RMS, indicates better tracking performance. The RMS 
of the control input R/S, was used as a measure for 
control effort. A higher RMS, indicates a higher control 
effort. 

A time-domain parameter estimation technique based 
on maximum likelihood estimation was applied in order to 
obtain the parameters of the human manual control model 
(Eq. (2)) [11]. Initial parameter estimates were obtained 
using a genetic algorithm, and refined using a gradient- 
based Gauss-Newton estimation. Manual control behav- 
ior was characterized by the visual error position gain K,,, 
lead time constant T7,, visual error velocity gain K,T7y,, 
time delay 7,, neuromuscular frequency w,,, and neuro- 
muscular damping ratio ¢,,. The gains kK, and K,T7, were 
of particular interest, since they indicate the relative use of 
error signal position and velocity information, which was 
expected to vary the most with different eccentricity an- 
gles. 
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The variance accounted for (V AF’) is a measure of 
how much of the measured control input signal u was ex- 
plained by the linear transfer function H,. It is an indi- 
cation of how linearly the human operator behaves, and is 
typically used as a measure for the goodness of fit of the 
linear human operator model [11]. 


Finally, the open-loop crossover frequency w,. and 
phase margin y,,, were determined as measures for track- 
ing performance and stability in the frequency domain. 


Hypotheses 


The main hypotheses of this research revolve around one 
question: what are the effects of viewing visual informa- 
tion from a PFD at different retinal eccentricities on hu- 
man manual control behavior? Since past research only 
investigated the effects of retinal eccentricity in tasks with 
human subjects as passive observers, the hypotheses for- 
mulated here cannot be based solely on findings from 
these previous studies. 


It was hypothesized that increasing the horizontal ec- 
centricity angle would make it harder for the participants 
to discern the error position. Therefore, it was expected 
that the visual error position gain K,, would decrease (hy- 
pothesis H1). In addition, previous research found veloc- 
ity perception thresholds increased with larger eccentric- 
ity angles. It was expected the visual error velocity gain 
k,Tr, would increase for increasing eccentricity angles 
(H2), as participants would need higher gain control in- 
puts to observe a change in the error. 


However, previous research also found that perceived 
velocity decreased with larger eccentricity angles, pos- 
sibly introducing a limit on the magnitude of K,T'y at 
higher eccentricity angles. Therefore, the expected in- 
crease in the error velocity gain KT, (as hypothesized 
in H2) was expected to be reduced or even canceled at 
higher eccentricity angles (H3). As participants were ex- 
pected to put a higher gain on error velocity information 
(higher K,T7,) for DI-like dynamics compared to SI-like 
dynamics, this effect was expected to be more dramatic 
for DI-like dynamics. 


Furthermore, conditions C2 and C6, and C4 and C8, 
had the same angle of eccentricity; however, in the hori- 
zontal or vertical axes, respectively. Since the task was a 
pitch control task, it would be harder to perceive the zero 
reference line of the fixed aircraft symbol, and thus the 
magnitude of the error e, when the PFD was below the 
fixation point. Therefore, it was hypothesized that perfor- 
mance would be lower for conditions C4 and C8 (H4). 
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Table 3 Main analysis of variance results. 


dependent DYN POS DYN x POS 
variable df F Sig. df F Sig. af F Sig. 
RMS, 1.0,9.0 55.594 xe 1.6, 14.699 46.707 ne 3.0, 27.0 3.773 ae 
RMS, 1.0, 9.0 2.711 - 3.0, 27.0 1.943 - 3.0, 27.0 2.289 - 
VAF 1.0, 9.0 13.156 ‘ch 3.0, 27.0 18.796 cl 3.0, 27.0 5.328 ie 
Kp 1.0,9.0 50.766 bes 3.0, 27.0 51.725 em 3.0, 27.0 21.183 ve 
TL 1.0,9.0 39.281 ae 3.0, 27.0 11.686 ee 3.0, 27.0 3.780 ad 
Kole 1.0, 9.0 9.273 ia 3.0, 27.0 1.958 - 1.5, 13.799 7.813 ** 
Ty 1.0, 9.0 0.829 - 3.0, 27.0 29.709 ¥e 3.0, 27.0 2.689 * 
Wn 1.0,9.0 21.840 n% 3.0, 27.0 5.171 se 3.0, 27.0 0.573 - 
(ee 1.0, 9.0 0.971 - 14,112,928 4.872 is L110074 4.324 * 
We 1.0,9.0 24.643 is 3.0, 27.0 50.577 *% 3.0, 27.0 8.085 oa 
Qm 1.0,9.0 69.875 *# 3.0, 27.0 13.369 #% 3.0, 27.0 3.335 ha 
** = — significant (p < 0.05) 
* = marginally significant (0.05 < p < 0.1) 
_ - not significant (p > 0.1) 
99 =  Greenhouse-Geisser sphericity correction 


RESULTS 


In this section, significant effects on the dependent vari- 
ables are discussed. The black solid lines depict data from 
the conditions with single-integrator-like dynamics (SI), 
and the grey dotted lines represent the conditions with 
double-integrator-like dynamics (DI). The error bars in- 
dicate the 95% confidence intervals of the mean for all 
participants, corrected for between-subject variability. 
Table 3 provides the results of the two-way repeated 
measures analysis of variance (ANOVA) performed on 
the dependent variables of the experiment. As part of 
the analysis, checks for outliers, normal distribution, and 
homogeneity of variances were performed. Normality of 
the data was assessed by Shapiro-Wilk’s test of normality. 
Homogeneity of variances was assessed by Levene’s test 
of equality of variances. As dependent variables contained 
no outliers and only some were non-normally distributed 
in only a few conditions, no corrections were applied to 
the data. Whenever the assumption of homogeneity of 
variances was violated, a Greenhouse-Geisser correction 
was used for the degrees of freedom of the F-distribution. 
For each dependent variable, statistically significant inter- 
actions between dynamics and display position are dis- 
cussed first, if they existed. Statistically significant main 
effects are discussed only if no statistically significant in- 
teractions were found. ANOVA results for simple main 
effects and post-hoc tests are given throughout the text. 


Tracking Performance and Control Activity 


Fig. 4a shows participants’ tracking performance in terms 
of the RMS of the error signal. A lower RMS, 
means a higher performance. There was a significant 
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two-way interaction between dynamics and position (Ta- 
ble 3). Performance was statistically significantly better 
for easier SI-like dynamics than DI-like dynamics for all 
four display positions (F'(1,9) = 36.123,p < 0.001, 
F(1,9) = 32.924,p < 0.001, F(1,9) = 33.984,p < 
0.001, F(1,9) = 28.937,p < 0.001). In addition, 
statistically significant differences in performance were 
found between display positions for both SI-like dynam- 
ics (F'(3, 27) = 47.665, p < 0.001), and DI-like dynam- 
ics (F'(3,27) = 21.848,p < 0.001). Post-hoc analysis 
with Bonferroni adjustment revealed that for the SI-like 
dynamics performance significantly degraded from foveal 
to +15 deg horizontal eccentricity (p < 0.001) and from 
+15 deg to +30 deg horizontal eccentricity (p = 0.005). 
Performance was equal between the +30 deg horizontal 
and -15 deg vertical eccentricity conditions (p = 1.000). 
For the DI-like dynamics, performance decreased from 
F to R (p = 0.005), but was equal between R and RR 
(p = 0.279) and RR and B (p = 1.000). 

The RMS of the control input is depicted in Fig. 4b. A 
higher RMS, means a higher control activity. No sta- 
tistically significant interaction or statistically significant 
main effects were found for any of the factors. However, 
Fig. 4b might suggest that control activity was slightly 
higher for the DI-like dynamics. 


Variance Accounted For 


The VAF is depicted in Fig. 5. A significant inter- 
action between DYN and POS was found. The VAF 
was equal between the two controlled dynamics for the 
foveal and the +15 deg horizontal eccentricity positions 
(F(1,9) = 3.533, p = 0.093 and F'(1,9) = 0.281, p = 
0.609, respectively). For the +30 deg horizontal eccen- 
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Figure 4 Tracking performance and control activity. 
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Figure 5 Variance accounted for. 


tricity position and the -15 deg vertical eccentricity po- 
sition, the VAF was statistically significantly higher for 
the DI-like dynamics (F'(1,9) = 10.959, p = 0.009 and 
F(1,9) = 34.382,p < 0.001, respectively), indicating 
that participants behaved more linearly than for the SI- 
like dynamics. For both the SI-like and the DI-like con- 
trolled dynamics, significant differences were found be- 
tween display positions (F'(3, 27) = 20.094,p < 0.001 
and F'(3, 27) = 8.439, p < 0.001, respectively). For the 
SI-like dynamics, the VAF was equal between the foveal 
and +15 deg horizontal eccentricity (p = 0.365), but then 
decreased between +15 and +30 deg horizontal eccentric- 
ity (p = 0.007). The VAF in the -15 deg vertical ec- 
centricity condition was equal to the VAF in the +30 deg 
horizontal eccentricity condition (p = 1.000). The VAF 
shows a similar trend for the DI-like dynamics. The VAF 
was equal between F and R (p = 0.113), was significantly 
lower in RR compared to F (p < 0.001), and significantly 
lower in B compared to F (p = 0.010). The VAF was 
equal between R and RR (p = 1.000), and RR and B 
(p = 1.000). 
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Manual Control Behavior 


Fig. 6a depicts the human operator gain on visual error 
position. A significant two-way interaction was found 
between DYN and POS (Table 3). The visual gain 
was Statistically significantly lower for the DI-like con- 
trolled dynamics for all display positions (F(1,9) = 
88.294,p < 0.001, F(1,9) = 35.207,p < 0.001, 
F(1,9) = 22.481,p = 0.001, F(1,9) = 14.700,p = 
0.004). Statistically significant differences were intro- 
duced between display positions for the SI-like dynamics 
(F(3, 27) = 49.979, p < 0.001), and the DI-like dynam- 
ics (F'(3,27) = 18.034,p < 0.001). Post-hoc analysis 
with Bonferroni adjustment revealed that for the SI-like 
dynamics the gain for position R is significantly smaller 
than for position F (p < 0.001), the gain in RR is smaller 
than in R (p = 0.011), and the gains in RR and B are 
equal (p = 1.000). For the DI-like dynamics, the gain 
was found to marginally decrease from F to R (p = 0.076) 
and from R to RR (p = 0.057), and then significantly in- 
creased from RR to B (p = 0.009). The visual gain was 
equal between R and B (p = 1.000). 

The visual lead time constant is depicted in Fig. 6b. A 
statistically significant two-way interaction was found be- 
tween controlled dynamics and display position. A very 
similar, but opposite trend was observed in T7, compared 
to K,. The visual lead time constant was higher for the 
DI-like dynamics for positions F, R, and RR (F'(1,9) = 
169.023,p < 0.001, F(1,9) = 35.469,p < 0.001, 
F(1,9) = 21.032,p = 0.001), indicating an increased 
reliance on visual rate information. For the -15 deg ver- 
tical eccentricity position (B), the visual lead time con- 
stant was equal between the two controlled dynamics 
(F(1,9) = 1.113,p = 0.319). For the SI-like dynam- 
ics, significant differences were found between different 
display positions (F'(1.4,12.9) = 7.553,p = 0.011). 
An increasing trend can be observed going from F to 
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Figure 6 Human controller parameters. 


B; however, J’; is found to be statistically significantly 
higher in RR compared to F and R only (p = 0.008 
and p = 0.017). Significant differences between dis- 
play positions were also found for the DI-like dynamics 
(F(3,27) = 8.147,p = 0.001). In this case, Ty, was 
found to be statistically significantly higher in RR com- 
pared to F and B only (p = 0.004 and p = 0.033). 


The overall human operator gain on visual error rate 
(K,T7) is depicted in Fig. 6c. Note that this gain is just 
the multiplication of kK, and Ty,. A significant two-way 
interaction was found for this gain as well. For the foveal 
and +15 deg horizontal eccentricity positions, the error 
rate gain was significantly higher for the DI-like dynamics 
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(F(1, 9) = 30.722, p < 0.001 and F'(1, 9) = 14.095, p = 
0.005). For the +30 deg horizontal and -15 deg verti- 
cal eccentricity positions, K,77, was equal between both 
controlled dynamics (F'(1,9) = 3.320,p = 0.102 and 
F(1,9) = 0.290,p = 0.603). Statistically significant 
differences were found between the different display po- 
sitions for both SI- and DI-like dynamics (F'(3,27) = 
3.505, p = 0.029 and F'(3,27) = 10.108,p < 0.001). 
Post-hoc analysis with Bonferroni adjustment indicated 
there were no statistically significant differences between 
any display position pairs for the SI-like controlled dy- 
namics. For the DI-like dynamics, post-hoc analysis re- 
vealed that 77, was significantly lower in R and B com- 
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pared to F (p = 0.008 and p = 0.004). K,7T7, was equal 
between F and RR (p = 0.194), and R and B (p = 1.000). 

There was no significant two-way interaction intro- 
duced in the visual time delay (Fig. 6d). However, a sta- 
tistically significant main effect introduced by the display 
position was found (Table 3). The visual time delay was 
equal for both dynamics and increased from the foveal 
to the +15 deg horizontal eccentricity display positions 
(p = 0.011), increased from the +15 deg to the +30 deg 
horizontal eccentricity positions (p = 0.008), and finally 
remained constant between the +30 deg horizontal and -15 
deg vertical eccentricity positions (p = 1.000). The dif- 
ference between the lowest time delay in F and the highest 
time delay in RR and B was approximately 100 ms. 

No statistically significant two-way interaction between 
DYN and POS was introduced in the human operator’s 
neuromuscular frequency (Fig. 6e). However, the main ef- 
fects of both DYN and POS were statistically significant. 
The neuromuscular frequency was significantly higher for 
the SI-like dynamics compared to the DI-like dynamics 
(Table 3). Post-hoc analysis with Bonferroni adjustment 
revealed that the neuromuscular frequency in RR was sta- 
tistically significantly lower than in F (p = 0.025), but 
was similar between R and F (p = 0.200), and B and 
F (» = 0.078), and more similar between R and RR 
(p = 1.000), and R and B (p = 1.000). 

Fig. 6f depicts the neuromuscular damping ratio. There 
was no significant two-way interaction introduced in this 
variable. The main effect of display position was statis- 
tically significant. Post-hoc analysis with Bonferroni ad- 
justment revealed that the damping ratio was statistically 
significantly higher in RR compared to F (p = 0.035), but 
was similar between R and F (p = 0.354), and B and F 
(p = 0.260). The neuromuscular damping ratio was also 
similar between RR and R (p = 0.279), and RR and B 
(p = 1.000). 


Open-Loop Characteristics 


Fig. 7 shows the open-loop crossover frequencies and 
phase margins. In general, when tracking performance 
increases, as indicated by an increase in crossover fre- 
quency, stability margins decrease; that is, performance 
is increased by sacrificing stability in the control loop. 
This effect was also observed in this experiment. A 
statistically significant two-way interaction between the 
controlled dynamics and the display position was found 
in the crossover frequency (Fig. 7a) and phase margin 
(Fig. 7b). The crossover frequency was statistically sig- 
nificantly lower for the SI-like dynamics for all display 
positions (F'(1,9) = 9.049,p = 0.015, F(1,9) = 
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14.331,p = 0.004, F(1,9) = 23.203,p = 0.001, 
F(1,9) = 41.458,p < 0.001). The phase margin was 
statistically significantly higher for the SI-like dynamics 
for all display positions (F'(1,9) = 153.406, p < 0.001, 
F(1,9) = 49.453,p < 0.001, F(1,9) = 14.800,p = 
0.004, F'(1, 9) = 31.037, p < 0.001). 

For the both the SI-like and DI-like dynamics, the 
crossover frequency was statistically significantly differ- 
ent between display positions (F'(3,27) = 52.209,p < 
0.001 and F'(3, 27) = 18.417, p < 0.001). For the SI-like 
dynamics, the crossover frequency significantly decreased 
from F to R (p < 0.001), and from R to RR (p = 0.012), 
and was equal between RR and B (p = 1.000). For the 
DI-like dynamics, the crossover frequency significantly 
decreased from F to R (p = 0.002), but remained ap- 
proximately constant between R and RR (p = 0.225), and 
RR and B (p = 0.247). The phase margin for the SI-like 
dynamics was equal between F and R (p = 0.872) in- 
creased from F to RR (p = 0.020), and was equal again 
between RR and B (p = 1.000). For the DI-like dy- 
namics, the phase margin was similar between F and R 
(p = 1.000), and F and B (p = 0.575); however was sig- 
nificantly higher in RR compared to F (p < 0.001). 


DISCUSSION 


Human manual control behavior in an active control task 
was successfully analyzed using a cybernetic approach in 
experimental conditions with visual stimuli presented at 
different retinal eccentricities and with two different con- 
trolled dynamics. Manual control behavior adapted signif- 
icantly as visual stimuli were presented at different eccen- 
tricity angles, as observed by significant variations in the 
parameters of the human operator model. Variations in the 
human control parameters were similar for both controlled 
dynamics, except for the parameters relating to visual rate 
information (TJ, and K,T7). 

The human operator visual gain kK’, decreased with in- 
creasing horizontal eccentricity for both controlled dy- 
namics, as was hypothesized in H1. This reduction in vi- 
sual gain was most likely caused by the the fact that the 
visual error was more difficult to discern as the horizon- 
tal eccentricity increased. The visual lead time constant 
Ty increased with increasing horizontal eccentricity, in- 
dicating that human controllers relied more on visual rate 
information as it became harder to observe the error. Com- 
bined, these results reflect the fact that velocity is gener- 
ally easier to perceive in peripheral vision compared to 
position. 

An interesting interaction was observed for the error 
velocity gain K,Ty, in Fig. 6c. For the SI-like dynam- 
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Figure 7 Open-loop performance and stability. 


ics, KT, slightly increased with increasing horizontal 
eccentricity, as hypothesized in H2. However, for the 
more difficult DI-like dynamics, K,7T7, decreased when 
using peripheral vision. Furthermore, KT, was equiv- 
alent between both controlled dynamics for the +30 deg 
horizontal and -15 deg vertical eccentricities. Minimum 
velocity perception thresholds are higher with increased 
retinal eccentricity [6]. As a result, and with the increased 
reliance on visual rate information, participants had to in- 
crease their gain on visual rate information in order to see 
changes in visual stimuli. This is reflected in the increase 
in KT, for the SI-like dynamics. However, at the same 
time, previous research also found that perceived veloc- 
ity decreased for increasing eccentricities [2]. It might be 
that this introduced a limit on the error velocity gain in 
peripheral vision, resulting in the reduction of K,Ty, for 
the DI-like dynamics (as hypothesized in H3); that is, as 
K,Ty, was already higher for the DI-like dynamics com- 
pared to SI-like dynamics for foveal vision and above the 
limit for peripheral vision, it could only go down to the 
limit value for increasing eccentricities. 


The visual time delay and neuromuscular frequency 
proved to be the human operator limitations most affected 
by retinal eccentricity. The visual time delay increased 
with increasing horizontal eccentricity, as also found in 
previous research [6]. The neuromuscular frequency de- 
creased with increasing horizontal eccentricity, indicating 
a reduction of the maximum frequency of visual changes 
acted upon by participants. 


The VAF was significantly lower for higher horizontal 
eccentricities and more so for the SI-like dynamics, indi- 
cating participants behaved more nonlinearly. There are a 
few possible explanations for this. First, higher frequen- 
cies are less visible at larger angles of eccentricity, there- 
fore the human controller will react less strongly to these 
frequencies (as observed in the neuromuscular frequency), 
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making his/her behavior more nonlinear. Second, since 
minimum velocity thresholds are higher at higher eccen- 
tricities, participants might have applied higher gain in- 
puts in order to perceive changes in error velocity. These 
additional inputs are not correlated with the target signal 
f; and will reflect as nonlinear behavior, decreasing the 
value of the VAF. 


The above-discussed changes in human manual con- 
trol parameters and linearity of control behavior when us- 
ing peripheral vision resulted in reduced tracking perfor- 
mance at higher eccentricities, as observed by higher val- 
ues for RMS, and lower values for w,. The reduction in 
performance was accompanied by an increase in stability 
margins. 


Conditions R and B had the same eccentricity angles, 
but in the horizontal and vertical axes, respectively. De- 
spite the fact that the eccentricity angles were equal be- 
tween these two conditions, the two conditions were quite 
different in nature because of the type of control task and 
the display used. First, for horizontal eccentricities, the 
zero-angle reference of the error was easier to discern than 
for vertical eccentricities, as it was in line with the fixation 
point. Second, for vertical eccentricities, the error moved 
towards and away from the fixation point, and, there- 
fore, the instantaneous eccentricity angle varied, whereas 
it was constant for horizontal eccentricities. Hypothesis 
H4 stated that the performance was expected to be worse 
in B compared to R, despite the equivalence in eccentric- 
ity angles. This was confirmed in Fig. 4a, where the worst 
performance was seen in RR and B for both controlled dy- 
namics. In fact, for most dependent variables B was more 
similar to RR, which indicates that in this pitch tracking 
task, the -15 deg vertical eccentricity condition had a big- 
ger impact on manual control behavior and performance 
than the +15 deg horizontal eccentricity condition. For 
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different control tasks and displays, these effects will most 
likely be different. 

From a practical standpoint, a few lessons can be 
learned from this study. The time delay, as observed in 
Fig. 6d, increased up to 100 ms with increasing eccentric- 
ities. This can be crucial during critical flight scenarios 
where reaction time is important, especially as reaction 
times in such scenarios can already be higher because pi- 
lots’ increased cognitive load. Second, the results showed 
that in pitch tracking tasks, viewing the PFD with a verti- 
cal eccentricity angle results in lower performance com- 
pared to viewing the task with the same horizontal ec- 
centricity angle. Furthermore, participants behaved more 
nonlinearly at higher eccentricity angles, indicating that 
pilots could possibly exacerbate unstable flight conditions, 
such as an aerodynamics stall, when viewing the PFD pe- 
ripherally. These results could be used to improve the de- 
sign of human-centered displays and cockpit layouts. 


CONCLUSIONS 


This paper investigated the effects of retinal eccentricity 
on human manual control behavior. Ten participants per- 
formed a pitch tracking task while looking at a simpli- 
fied primary flight display at different angles of horizon- 
tal and vertical eccentricity, and with two different con- 
trolled dynamics. Error position gain generally decreased 
with increased eccentricity. Error velocity gain increased 
with eccentricity for the single integrator-like dynamics 
and decreased for the double integrator-like dynamics, ef- 
fects which might be linked to higher minimum velocity 
perception thresholds and decreased perceived velocity in 
peripheral vision. Vertical eccentricity resulted in a larger 
impact on control behavior and lower performance than 
equivalent horizontal eccentricity, as control behavior and 
performance in the -15 deg vertical eccentricity condition 
was most similar to the +30 deg horizontal eccentricity 
condition. The human operator time delay was up to 100 
ms higher at higher eccentricities compared to foveal vi- 
sion. Furthermore, participants behaved more nonlinearly 
at higher eccentricity angles. The experiment provided 
important insights into human manual control and perfor- 
mance using peripheral vision, which could be considered 
in the design of human-centered displays and cockpit lay- 
outs. 
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